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Yo.5Cao.5BaCc-407 was recently introduced as a possible candidate for capturing some of the pre- 
dicted classical spin kagome ground state features. Stimulated by this conjecture we have taken up 
a more complete study of the spin correlations in this compound with neutron scattering methods 
on a powder sample characterized with high-resolution neutron diffraction and the temperature 
dependence of magnetic susceptibility and specific heat. We have found that the frustrated near 
neighbor magnetic correlations involve not only the kagome planes but concern the full Co sub- 
lattice as evidenced by the analysis of the wavevector-dependence of the short range order. We 
conclude from our results that the magnetic moments are located on the Co sub-lattice as a whole 
and that correlations extend beyond the two-dimensional kagome planes. We identify intriguing dy- 
namical properties, observing high-frequency fluctuations with a lorentzian line-width T < 20 meV 
at ambient temperature. On cooling a low-frequency (~ 1 meV) dynamical component develops 
alongside the high-frequency fluctuations, which eventually becomes static at temperatures below 
T w 50 K. The high-frequency response with an overall line-width of ~ 10 meV prevails at T < 2 K 
coincident with a fully elastic short-range ordered contribution. 

PACS numbers: 25.40.Fq, 61.05.Fm, 75.25+z, 75.40.Gb 
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I. INTRODUCTION 



An ensemble of magnetic moments ("spins") residing on 
a two-dimensional kagome lattice has been known for a 
long time as a prototype for topological magnetic frustra- 
tion!^] In recent years, a number of materials has been 
synthesized where 5=1/2 or 5=3/2 spins reside on such 
a lattice or have a sub-lattice containing kagome lay- 
ers. If the couplings are sufficiently two-dimensional in 
character, long-range order at low temperature is sup- 
pressed and the effects of frustration can be studied 
experimentally. Topical materials in this class include 
5=1/2 systems such as volborthite, Cu3V207(OH) 2 • 
2H 2 0,^I 5 and herbertsmithite, ZnCu 3 (OH) 6 Cl2P™ as 
well others with higher spin quantum number, such 
as the 5=5/2 jarosites AFe 3 (OH) 6 (S0 4 )2, liimj and the 
double kagome layer (termed pyrochlore slab) com- 
pounds with 5=3/2 - SrCr 9 ^Gai2-9*Oi9 [SCGO^PHIZ] 
and Ba 2 Sn 2 ZnCr 7;r Gaio-7x022 (BSZCGO).™^ Due to 
complicated chemical composition and crystal structure 
the low temperature properties of these systems are com- 
plex and sometimes controversial. Deviations from the 
simplest nearest neighbor isotropic interaction picture, 



slight deviations from the ideal frustrated topology of 
the magnetic sub-lattice and the role of defects and im- 
purities can influence the properties in a way that may 
be difficult to clarify. A long term goal in the field has 
been to synthesize compounds that would exhibit and 
therefore confirm the expected quantum (Ref. [22] and 
references therein) or classicaP^H] properties revealed in 
theoretical and numerical studies of the nearest neighbor 
Heisenberg hamiltonian. 

Recent studies of "114" cobaltite compounds with 
the composition YBaCo^^ (7 < x < 8) and sub- 
stituted variants have revealed the interesting frus- 
trated properties and intricate chemistry of these com- 
poundsJ^HSD The substitutional^ disordered system, 
Yo.5Cao.5BaCo407, 39 was proposed as a new kagome 
material involving Co 2+ (5=3/2) ionsP^ A combina- 
tion of different arguments involving charge balance, 
local Co ion site symmetry and possible oxygen defi- 
ciency was used to propound Yo.sCao.sBaCo^r as a 
two-dimensional spin system in which high-spin (5=3/2) 
Co 2+ ions are arranged in kagome planes (labelled Co2 
on the 6c Wyckoff position), with additional low-spin 
(5=0) Co 3+ ions preferentially residing in intermediate 
layers (labelled Col on the triangular planar 2a Wyck- 



off position), separating the kagome planes. Quantita- 
tive modeling of the reported cold neutron polarized neu- 
tron diffraction data led Schweika and co-workers to the 
conclusion that the observed features match the classi- 
cal kagome ordered state expectations^ allowing for the 
"weathervane" defects typical of the low-energy response 
of the system. 12 Motivated by this observation we have 
undertaken a more complete neutron scattering exami- 
nation aiming at better understanding of the static and 
dynamic magnetic response, in the context of theoreti- 
cal work on the dynamics and ordering of the classical 
kagome system ! 25 " 27 ! 



It was reported earlier that Yo.sCao.sBaC^Oy is char- 
acterized by a large magnitude Curie- Weiss constant, 
Ocw ~ —2200 K, and the apparent abse nce of a mag- 
netic ordering transition down to T ~ 1.2 Kp™2l suggest- 
ing a significant degree of frustration. Accordingly the 
magnetic properties of Yo.sCao.sBaCo^r appear very 
different from the parent compound YBaCo407 which 
shows a transition to a state with long range magnetic 
order below T = 110K, at a reasonably high tempera- 
ture, but nevertheless below the mean field expectation 
based on the value 0c w ~ —510 K found for this com- 
poundPSOne might argue that the difference with respect 
to the non-substituted case is associated with the im- 
portant effect of heterovalent substitution of Ca 2+ for 
Y 3+ and the resulting enhancement of geometric frus- 
tration due to decrease of dimensionality. The magni- 
tude of the Curie- Weiss constant indicates the expected 
energy scale of the dynamic response, and accordingly 
we have chosen to investigate the dynamic response of 
Y .5Cao.5BaCo407 over a broader energy range beyond 
the cold neutron regime investigated earlier. 40 



The presentation of our results is organized as follows. 
After a short account of the experimental techniques in 
section |ll] we present an overview of the characteristics of 
the sample examined in section III A We present a struc- 
tural study of Yo.5Cao.5BaCo407 using high-resolution 
neutron powder diffraction in section |IIIB[ then in the 
section [lire] we present new polarized neutron diffraction 
data, found to be in quantitative agreement with those 
reported earlier™ However, our analysis of the spatial 
spin correlations suggests that all Co ions carry moments 
and participate in the magnetic response. Furthermore, 



in section |IIID| we present a neutron spectroscopic inves- 
tigation which indicates a dynamical magnetic response 
reaching over a wide range of time and energy scales, and 
also revealing a freezing of the system when it is cooled 
to low temperatures T < 50 K, in agreement with our 
bulk magnetic observations. We show that due to the 
broad extent of the dynamics, polarized neutron diffrac- 
tion using cold neutrons misses a significant portion of 
the fluctuating moment, reinforcing the conclusion that 
all Co atoms are magnetic and that Yo.5Cao.5BaCo407 
is not a classical kagome system. 



II. EXPERIMENTAL METHODS 

A powder sample of Yo.5Cao.5BaCo407 (mass ^ 8 g) 
was prepared in a manner similar to that used previ- 
ously. 39 High-purity Y 2 3 , CaC0 3 , BaC0 3 , and Co 3 4 
were mixed in stoichiometric ratios and heated in air at 
1100 C with several intermediate grindings. After the 
final firing at 1100 C, the sample was quenched to room 
temperature to minimize the possibility of excess oxy- 
gen stoichiometry in the sample. This quenching from 
high temperatures should also lead to a random substi- 
tution of Ca for Y. X-ray powder diffraction was per- 
formed with Cu-K a radiation, showing that the sample 
was well crystallized and single-phase within the sensi- 
tivity of the technique. The refined cell parameters at 
room temperature were consistent with values reported 
earlier by Valldor, 39 although we shall see below that 
their reported structure (hexagonal P63171C space group) 
does not match our observations. 

Measurements of DC and AC susceptibility were per- 
formed in a Quantum Design MPMS between T = 2 K 
and 400 K. The specific heat data were taken in two sets, 
one over a range 4K<T<100K without applied field, 
another over a range 20 K < T < 300 K with applied 
field up to jaoH = 14 T. The lower T data were obtained 
using a semi-adiabatic method (MagLabVSM, Oxford In- 
struments) and were corrected for the sample holder and 
grease contributions. The high temperature specific heat 
from 20 K to 300 K was measured with an high precision 
continuous heating adiabatic calorimeter. 

The neutron scattering measurements were made at 
the Institut Laue-Langevin (ILL) in Grenoble, France, 
and the ISIS Pulsed Neutron facility in Oxfordshire, 
UK. A high-resolution neutron powder diffraction struc- 
tural study was performed at the HRPD diffractometer 
at ISIS. An xyz neutron polarization analysis diffrac- 
tion study was performed using the D7 diffuse scatter- 
ing instrument at the ILL, with an incident wavelength 
of 4.8 A. Time-of-flight (TOF) inelastic neutron spec- 
tra were measured using the thermal spectrometer IN4 
at the ILL, with various settings of the incident wave- 
length (Xi = 1.1 A, 1.3 A, 2.2 A) and instrumental reso- 
lution in the range between A(W) « 3.5 meV to 0.7 meV 
(FWHM). Additional inelastic TOF data with higher res- 
olution A(Huj) « 100 /xeV (FWHM) were taken using 
the IN5 cold neutron spectrometer at the ILL operat- 
ing at Xi = 5.0 A. Elastic scans with 1 fieV resolution 
were recorded on the backscattering instrument IN 16 at 
the ILL. In these experiments, data were corrected for 
neutron detection efficiency and background scattering 
using standard procedures. For further investigation of 
the slow dynamics a neutron spin echo (NSE) experi- 
ment was performed using the multidetector spin echo 
spectrometer IN11C at the ILL, with a neutron wave- 
length of 5.5 A (bandwidth 15%), giving a dynamic range 
from 5 x 10 -12 s to 2 x 10 ~ 9 s and simultaneous cover- 
age of a Q range from 1 A -1 to 1.6 A -1 . A sample of 
Hoo.7Yi. 3 Ti2 07 measured at 1.5 K was used to measure 
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FIG. 1: (Color online) The temperature dependent DC mag- 
netic susceptibility at fioH = 100 G (main figure). Zero-field- 
cooled (ZFC) and field-cooled (FC) branches are shown The 
inset shows the shift to higher temperature of the cusp in the 
AC susceptibility as a function of frequency, measured using 
a primary AC field of /noH = 3 G (rms) for / = 0.1 (black), 
1 (blue), 10 (red) and 100 Hz (green). 



the instrumental resolution in this experiment. 41 

D7 and IN11C use cold neutrons with low incident en- 
ergy {E{ ~ 3.5 meV in the experiments reported here) 
and the supermirror analyzer devices used on these spec- 
trometers effectively limit the energy transfer acceptance 
for neutron energy gainP^ This means that the effective 
energy transfer window of these measurements covers the 
range -10 meV < ftw < 3 meV. 



III. EXPERIMENTAL RESULTS 

A. Susceptibility and specific heat 

Fig. [T] shows the temperature dependent DC magnetic 
susceptibility obtained at jiqH = 100 G. The main fea- 
ture is a sharp cusp in the zero field cooled (ZFC) sus- 
ceptibility at a temperature of T g = 48 K concomitant 
with a bifurcation of the ZFC and field cooled (FC) sus- 
ceptibilities, indicative of spin freezing. 43 At the same 
temperature the AC susceptibility shows a characteris- 
tic, weakly frequency dependent cusp as often observed in 
systems of freezing spins. There are similarities with the 
data reported earlier by Valldor, 39 but our susceptibil- 
ity data has fewer peaks and is much less structured and 
smoother (there were at least three clear peaks in Vall- 
dor 's data set). The third anomaly reported by Valldor 
at T = 14 K is missing, indicating that it may depend on 
the preparation of the specific Yo.sCao.sBaCo^r sam- 
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FIG. 2: (Color online) The temperature dependence of the 
specific heat, (top panel) Data measured in zero magnetic 
field of between 20 K < T < 300 K (blue, vertical lines), and 
in a magnetic field of fi H = 14 T between 20 K < T < 120 K 
(green, crosses). Note that these two are totally overlapping 
indicating no field dependence over the T-range concerned. 
The red points show the low-T portion in zero field, 
(lower panel) The low-T data in zero field (blue solid circles) 
fitted to a Debye phonon-term and an exponentially activated 
contribution oc exp(— i^/Zc^T) (blue solid line) as well as the 
extracted purely Debye part (blue dashed line) and the differ- 
ence between the data and the Debye contribution (red hollow 
circles) , with the exponential part of the fit (red solid line) . 



pie. There is no evidence of Curie- Weiss like behavior 
up to the maximum temperature measured. 

Fig. [2] shows the measurement of the specific heat at 
zero field and an applied field of fioH = 14 T. No visible 
anomaly appears in the vicinity of T g with or without 
an applied field, suggesting the absence of a true col- 
lective thermodynamic transition. Hence, no long range 
order or any low-energy magnon-like excitations asso- 
ciated with such a state can be expected at low tem- 
perature T < T g . The trend as T — »> 300 K suggests 
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FIG. 3: (Color online) Neutron powder diffraction pattern of 
Yo.5Cao.5BaCo407 at room temperature, measured using the 
HRPD diffractometer at ISIS. A Rietveld refinement model 
of the data, calculated using FULLPROF 44 ^ is shown as a solid 
line through the data. The difference between the data and 
model is shown as a solid line below. 



that the specific heat overshoots the high-temperature 
Dulong-Petit limit (3 &b per ion) of the lattice specific 
heat which coincides with the top of the vertical scale 
of the plot (324 J mol _1 K _1 ). The specific heat at low 
temperature, T < 20 K, is best fitted using two terms; 
a main contribution that follows a Debye T 3 dependence 
and an activated term oc exp(— Ea/ksT). Fitting the 
low-temperature data to this model we find a Debye 
temperature Od ~ 290 K and an activation energy of 
^aA#=25...30 K, depending on the fitting range. As a 
hypothesis we can consider that the exponential increase 
of the specific heat at low temperature and the overshoot 
of the Dulong-Petit limit at high temperature are due to 
a magnetic contribution to the specific heat. We shall re- 
turn to this issue below with regard to the the dynamic 
response obtained from our neutron spectroscopic data. 



B. High-resolution neutron powder diffraction 
study 



TABLE I: Atomic parameters for Yo.sCao.sBaC^O? at room 
temperature. The symmetry is trigonal P31c, with cell pa- 
rameters, a = 6.3075(1) A and c = 10.2388(1) A. The r- 
factors were found to be R w = 6.05, i?Bra gg = 5.52 and a 
goodness of fit (GoF) = R w /R exp = 1.83. 
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FIG. 4: (Color online) Crystal structure of Yo.5Cao.5BaCo407 
at room temperature. Y/Ca atoms in purple, Ba in grey, Col 
centered tetrahedra in blue, Co2 centered tetrahedra in red, 
01 in yellow/green, 02 in blue and 03 in red/pink. The 
splitting of the 01 and 03 oxygen sites along the Ba-0 bonds 
causes the base of oxygen tetrahedron around the Co2 posi- 
tion (red) to be either parallel to the kagome plane, or to be 
tilted by an angle of 8.7°. 



A high-resolution neutron powder diffraction pattern 
of Yo.5Cao.5BaCo407 is shown in Fig. [3] Rietveld analy- 
sis of the data was carried out using the FULLPROF pro- 
gram™ The best model fit to the data was found for 
the trigonal P31c space group, in agreement with most 
other studies of "114" series cobaltites at high tempera- 
^ ures |33|34J45j fo u i m contrast to the PQsmc space group 
found by Valldor using Cu-K a x-ray radiationPS This 
is likely due to the fact that x-ray diffraction - which is 
much less sensitive to oxygen than neutron diffraction - 
was used to solve the structure in that case. 



The cell and atomic parameters found from the Ri- 
etveld refinement are given in Table [T| and the crystal 
structure is depicted in Fig. [4] Evident from the ta- 
ble is that there are 5 distinct oxygen sites with rea- 
sonable thermal parameters in the unit cell, as opposed 
to 3 in the case of pure YBaCo407 with a splitting of 
the 01 (6c) and 03 (6c) sites along the Ba-0 bonds. For 
Yo.5Cao.5BaCo407 large anomalous thermal parameters 
were found for oxygen atoms in the P63171C model,^ in- 
dicating some degree of oxygen disorder. The occupancy 



of the split 01 and 03 positions was found to be ap- 
proximately in proportion to the Y-Ca occupation ratio, 
with the refined Y-Ca occupation found to be 0.55 for 
Y and 0.45 for Ca. Thus we surmise that the choice of 
01 and 03 position in the unit cell depends on whether 
that position is coordinated with a Y or Ca ion. The 
Y-Ca disorder in the lattice could well have such an ef- 
fect on the crystal structure, given that; a) the ionic 
radii of Y 3+ (1.04 A) and Ca 2+ (1.14 A) in an octa- 
hedral site are different by ~ 10 %, and; b) in the pure 
Y-114 sample the C0O4 tetrahedra are aligned to the 
kagome planes J^ while in the pure Ca-114 sample the 
C0O4 tetrahedra form a corrugated structure out of the 
kagome plane. 46 In Yo.sCao.sBaCo^r the C0O4 tetra- 
hedron around the kagome Co2 site is either orientated 
with basal plane aligned with the kagome plane (Oil and 
031 positions) - as in the pure Y sample, or oriented with 
the basal plane at an angle of 8.7° to the kagome plane 
(012 and 032 positions) - as in the pure Ca sample. 
Given that the ratio of site occupation of the 011:012 
and 031:032 sites is roughly in proportion to the Y:Ca 
occupation ratio, we conclude that the C0O4 tetrahedra 
are tilted when the apex is connected to a Ca ion, and 
aligned when connected to a Y ion. The Co-0 distances 
fall in the range 1.84-1.99 A and 1.89-2.01 A for the Col 
site and Co2 site C0O4 tetrahedra respectively. These 
distances compare well with those found in YBaCo407 
and CaBaCo 4 7 . [1SMI Similarly to YBaCo 4 7 , the C0O4 
tetrahedra at the Co2 site are larger than those at the 
Col site, suggesting that the kagome layer of tetrahedra 
should be preferentially Co 2+ . Another striking feature 
of the P31c model concerns the Ba coordination. The 
twelve nearest oxygen neighbors are located at distances 
typical for barium, falling in the range 2.71 A to 2.99 A. 
In other 114 cobaltites, most of the Ba-0 distances are 
more th an 3 A and cannot be classified as true Ba- 
O bonds )2£l 45 l 46 l In these cases, Ba 2+ is severely under- 
bonded and it has been proposed that some 114 cobaltites 
undergo a trigonal-orthorhombic structural transition to 
compensate for this underbondingP^ With regard to oxy- 
gen stoichiometry, the current data do not allow an accu- 
rate estimation. However, we can say that allowing the 
oxygen stoichiometry to be a free parameter in the Ri- 
etveld refinement did not improve the fit. Therefore our 
neutron diffraction data are consistent with an oxygen 
stoichiometry of 7. 



Concerning the Co2 kagome plane geometry, it is im- 
portant to note that the geometry is not perfect kagome, 
but rather a lattice of two distinct sizes of corner-sharing 
equilateral triangle. Indeed, a perfect kagome lattice is 
only possible for a general position of type (x, 1 — x, y) 
where x = 1/6. The side lengths of these two equilateral 
triangles making up the kagome planes were found to be 
3.141 A and 3.168 A. The Col-Co2 intersite distance is 
bigger (3.23 A). The distance between the kagome planes 
is found to be 5.12 A. 
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FIG. 5: The temperature and wavevector transfer dependence 
of the magnetic scattering (separated from nuclear contribu- 
tions using xyz polarization analysis) of Yo.sCao.sBaCo^r 
between room temperature and 3 K. 



C. xyz Neutron polarization analysis study 

Using the technique of xyz neutron polarization anal- 
ysis,^ which allows the separation of the magnetic neu- 
tron scattering cross-section from nuclear and nuclear 
spin contributions, we have measured the magnetic neu- 
tron scattering from Yo.sCao.sBaQ^Or at temperatures 
between 300 K and 3 K. Fig. [5] shows the evolution of 
the magnetic scattering as a function of temperature and 
neutron wavevector transfer, Q = k\ — k$ where k\ and 
kf are the incident and final neutron wavevectors respec- 
tively. 

The main feature in neutron diffraction is the appear- 
ance of a broad peak of diffuse scattering centered around 
\Q\ ~ 1.4 A -1 , indicating short ranged spin correlations. 
These correlations begin to appear at rather high tem- 
perature, with evidence of short-range order at the high- 
est temperature measured (300 K). The magnetic inten- 
sity appears to increase monotonically as the tempera- 
ture decreases. This indicates that the Co moments are 
becoming increasingly "static" on the timescale of the D7 
measurement. As explained in section [IT] the neutron 
cross-section measured on D7 is integrated over energy- 
transfer between the limits -10 meV and +3 meV, im- 
plying that the spin-fluctuation rate of the Co moments 
must be less than ~ 10 12 Hz if they are to contribute to 
the measured scattering cross-section on D7. The inte- 
grated magnetic moment derived from the temperature 
scan on D7 is shown in the inset of Fig. [6] Sharp Bragg- 
peaks may also be seen in the data at |Q| = 1.25 A -1 and 
\Q\ = 2.45 A -1 . These peaks, which are independent of 
temperature, are the first two magnetic Bragg peaks of 
CoO, which is an antiferromagnet below Tn = 290 K. 
The intensity of the CoO magnetic peaks, in conjunction 
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FIG. 6: The upper panel shows the experimentally deter- 
mined scattering cross section S(Q) at T=3 K together with 
an RMC model fit (solid line) assuming non-magnetic Col 
atoms. The inset shows the temperature dependence of the 
Co moment per formula unit as determined by our RMC cal- 
culations. The lower panel shows the real space, shell to shell 
spin correlations over distance. The solid line is a guide to 
the eye. 



with the fact that CoO peaks are absent from the pow- 



der neutron diffraction data presented in section IIIB[ 
allows us to put an upper limit of 1% by weight for the 
CoO impurity level. 

Modeling of the data was undertaken by performing a 
reverse Monte Carlo (RMC) fit of the D7 data. Initially, 
assuming a non-magnetic Col (2a) site, as was proposed 
in the previous study, these moments were excluded from 
the simulation. A box of 8 x 4 x 4 unit cells (making a 
total of 1536 Co2 atoms) was set up with an initially 
random spin configuration. In the simulation, the spins 
were assumed to be fully classical (Heisenberg 3D) and of 
uniform magnitude. In each Monte Carlo step, one spin 
was randomly chosen and rotated by a random angle, 
and the powder average of the cross section (dcr/dQ)M&g 
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FIG. 7: The correlation distributions in the first four near 
neighbor shells from RMC simulation at T = 3 K, assum- 
ing non-magnetic Col-atoms. The first two shells are in the 
kagome plane and show a strong 120° correlation 



was computed using 
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where the magnetic pre-factor 2/3(7 n rg/2) 2 = 0.049 
barns, f 2 (Q) is the Co (2+) squared magnetic form fac- 
tor, ji 2 is the squared magnetic moment per Co atom, r n 
and Z n are the distances and coordination numbers of the 



near neighbor Co coordination shell and ( Sq • S n 



,th 



is the average Heisenberg spin-correlation between unit 
spins on a central Co atom at the origin and the n th 
near neighbor Co coordination shell. This model was 
then compared with the data, and a x 2 goodness of fit 
parameter calculated. Each RMC step (individual spin 
rotation) was accepted if it resulted in an improved fit, 
and rejected if not. In order to verify that the calcu- 
lated model spin configuration was independent of the 
initial random starting position, many individual RMC 
calculations were performed, and the results averaged. 

The scattering cross-section da/dfl resulting from the 
RMC is shown in the upper panel of Fig. [6] as a solid 
line. The goodness of fit parameter obtained for this 
model was found to be x 2 — 2.61. 

Fig. [7] shows the correlation distributions in the first 
four near neighbor shells resulting from the RMC cal- 
culation. The first two shells (Fig. [7| top left) show 
a peak in the correlation distribution function at —0.5, 
that is, an average angle of 120° between spins located 
in the same triangle. From the oscillation period of the 
radial correlations (lower panel of Fig. |6| an antiferro- 
magnetic wavevector of 1.22 A -1 is estimated, half way 
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FIG. 8: The upper panel shows the experimentally deter- 
mined scattering cross section S(Q) at T=3 K together with 
the RMC simulation assuming both Col and Co2 atoms are 
magnetic (solid line), compared with the previous fit from 
Fig. [6] with the Col site excluded (dashed line). The lower 
panel shows the real space, shell to shell spin correlations over 
distance. The solid line is a guide to the eye. 



between Q = 1.15 A -1 and Q^ = 1.33 A -1 (in the 
notation of Ref. [40j), which denotes an ordered antifer- 
romagnetic structure with uniform chirality, where the 
magnetic and atomic unit cells coincide (AF ordering 
wavevector &af = 0), and a staggered chiral structure, 
where the magnetic unit cell is expanded by y/3 x y/3 
with respect to the atomic unit cell, (&af(1/3, 1/3)). 

The use of xyz polarization analysis enables the scat- 
tering cross section to be easily converted to absolute 
units, allowing a determination of the average magnetic 
moment, by allowing /i 2 in Eq. [I] to vary. In this way, 
we obtain a squared moment of/i 2 = 21.6 ± 0.8 n 2 B per 
formula unit at T = 3 K. Note that this moment varies 
strongly with temperature, as shown in the inset of Fig. [6] 
This is a similar number to that found by Schweika and 
co-workers™ and is the basis of their conjecture that the 
Col (2a) sublattice has zero moment. Assuming this, 
then the remaining three Co2 (6c) atoms have a squared- 
moment of 7.2 ± 0.3 /ig each. If we assume that the 
Co2 atoms are frozen, and that there is no orbital con- 



tribution to the moment, then the measured spin per Co 
would be S = fi/g s = 1.34 ± 0.03. This is around 90 
% of the expected S — 3/2 (on the Co2 site) and there- 
fore accounts for the majority of the Co moment. If, on 
the other hand, we were to assume that a moment exists 
on the Col sublattice (and that the measured total mo- 
ment per formula unit were shared out evenly between 
the Co positions), we would find a squared moment of 
only 5.4 ± 0.2 /i| per Co atom, corresponding to a spin 
S = 1.16 ± 0.02 or only 64% of the expected S = 1.81 
(made up of 2.5 5=3/2 Co 2+ and 1.5 5=2 Co 3+ atoms 
per formula unit 39 ). This fraction of ~ 36% of miss- 
ing Co moment is, however, readily explainable without 
the necessity of claiming a low-spin 5=0 state on the 
Col sublattice. Any Co moment with a fluctuation rate 
> 10 12 Hz will not contribute to the measured scattering 
cross-section due to the limited energy integration win- 
dow of the D7 spectrometer, as discussed above. There- 
fore, the missing moment fraction could well be indica- 
tive of high-frequency magnetic fluctuations, outside the 
measurement window. There are many cases of residual 
high-frequency quantum zero-point fluctuations at low 
temperatures in frustrated and disordered magnets, and 
these high-frequency fluctuations often lead to a lower- 
than-expected integrated magnetic moment, when mea- 
sured using cold neutrons ! 10 * 48 * We shall see below that 
such high frequency fluctuations are indeed observed in 
Y .5Cao.5BaCo407. 

Given that there is no a priori justification for the as- 
sumpti on o f a low-spin 5=0 moment on the Col site (see 
Section IV), we now revisit the analysis of the magnetic 



diffuse scattering, and include these sites in the RMC 
model. It is now assumed that both Co sites have the 
same moment. The fitted magnetic structure factor, with 
the Col site included, is shown in Fig. [8| together with 
the previous model with the Col site excluded. The in- 
clusion of the Col site results in a fit of improved qual- 
ity, with x 2 — 1-69, indicating that our powder neu- 
tron diffraction data seems to indicate the presence of 
a Col moment. We conclude therefore that the observed 
magnetic scattering cross-section is not consistent with a 
planar kagome-like magnetic ground state, but rather a 
3d antiferromagnetically correlated structure, with mag- 
netic moments on all the Co sites. This is consistent 
with the observation of magnetic moments on all of the 
Co sites in ordered stoichiometric YBaCo407J 34 * 38 * and 
disordered off-stoichiometric YBaCc^Or.i!™ 



D. Dynamic magnetic response 

The high energy dynamic response of the system as 
probed by thermal neutron TOF spectroscopy on the 
IN4 spectrometer shows many features not reported ear- 
lier. An intense ridge of quasi-elastic magnetic scattering 
centered around Q = 1.4 A -1 , at the same position as 
the peak in the "static" diffuse scattering measured in 
the more restrictive dynamic range of D7, is seen in the 
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FIG. 9: (Color online) Intensity maps of the scattering func- 
tion S(Q,lj) as measured using IN4 with an incident neutron 
energy of Ei=17 meV at three different temperatures, T= 20, 
48 and 100 K, from left to right. The curved feature ema- 
nating from Q w 0.7 A -1 is due to the gap in the detector 
assembly at about 20 «10°. 



(Q, uj) plane as depicted in Fig. [9] This ridge extends over 
a broad energy range with a non-trivial temperature de- 
pendence which will be detailed below. The broad energy 
range over which the ridge is observed, and its diffuse na- 
ture, suggest that it is due to an ensemble of correlated 
magnetic "clusters" with a large size distribution and a 
macroscopic number of interacting degrees of freedom. 

We can examine the line-shape and temperature de- 
pendence of this magnetic scattering independently of 
any particular model by considering the following gen- 
eral expression for the scattering function, 



S(Q,u,T) = (1 
= (1 



' kT 



^)-Y'(Q )W) T) (2) 



Using this relationship the product of the wave-vector 
transfer dependent susceptibility and the spectral weight 
function, xq(T)Fq(uj,T), can be evaluated from the ex- 
perimental data. The spectral weight function is normal- 
ized, f F(u,T)du = 1, such that the temperature depen- 
dence of the integrated inelastic intensity is given by the 
Q-dependent susceptibility, xq{T)- Fig. 10 shows the 



temperature dependence of xq(T)Fq(uj,T) taken from a 
constant Q-cut through S{Q,uS) centered at 1.4 A -1 , on 
cooling from T=280 K to T=1.5 K. The continuous lines 
on the plot follow a single lorentzian line-shape of the 
form 



Fq(u,T) = 



1 



1 



n w 2 + Tj 



(3) 



which describes the data at high temperatures rather ac- 
curately. Furthermore we note that the observed nar- 
rowing of the line-shape as the temperature is reduced, 
appears to take place in the energy transfer range \Huj\ < 
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FIG. 10: (Color online) Temperature dependence of the ex- 
perimentally evaluated spectral response at Q — 1.4 A -1 de- 
fined in Eq. 15] Data from bottom up taken at 250 K (red), 
200 K (green), 150 K (yellow), 100 K (black), 70 K (pink), 
48 K (purple), 20 K (light blue) and 1.5 K (dark blue). The 
ranges of the data are limited on the energy gain side (negative 
energy transfer) by the thermal population/detailed balance 
factor and at T=1.5 K only the energy loss data is shown. 
Various instrument settings were used and combined. The 
shaded area covers the instrumental resolution range ± 3 meV 
for T > 100 K and ± 0.7 meV for T < 100 K. The solid lines 
show examples of fits to the data of a single broad lorentzian 
component at T = 250 K (red solid line) and T = 100 K (black 
solid line). At T = 100 K we note the onset of a narrow com- 
ponent that is examined in detail in Figs. |12| and |13| while the 
broad component stays virtually independent of temperature 
below 100 K. 
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FIG. 11: The temperature dependence of the line-width F 
(red hollow circles) and of the intensity x (blue circles) at a 
wavevector transfer of Q = 1.4 A -1 (the peak of the suscep- 
tibility), of the broad lorentzian component of the inelastic 
response measured on IN4. 



5 meV. Therefore we proceed by restricting our atten- 
tion to the high-frequency fluctuations, fitting only the 
broad part of the response outside of this range to obtain 
temperature dependent values of the broad lorentzian 
line-width, Tq and intensity \Q- These are depicted in 
Fig." 
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The line-width r Q at Q = 1.4 A" 1 falls from 
20meV at 300 K to - 10 meV at 100 K, and thereafter 
remains roughly constant with temperature, falling only 
slightly to ~ 7 meV at 1.5 K. The peak susceptibility, xq 
at Q = 1.4 A -1 associated with this broad component, 
increases slowly between 300 K and 1.5 K. 

In our judgement, the line-shape of the total response 
is complex and the existence of a distinct broad com- 
ponent is not justified by any self-evident arguments. 
Nevertheless, considering that both Tq and \Q are on ly 
weakly temperature dependent for T < 100 K, and that 
no great deviation in these parameters is observed around 
the spin freezing transition seen in our susceptibility data 
(Fig.Jl]), we proceed further to characterize the narrower 
central part of the response. In order to accomplish this, 
we use data taken over a smaller dynamic range, but with 
better resolution using the IN5 spectrometer, shown in 
Fig. [l2j Keeping the parameters of the broad compo- 
nent fixed (r = 10 meV and \Q ~ 1) at temperatures 
below 80 K, we can fit the narrow response with another 
lorentzian form using Eqs. [2] and [3J shown in Fig. [l2j 
The elastic intensity and the integrated low-frequency 
response, found from the fitted IN5 spectra as a func- 
tion of Q and T are shown in Fig. [13] "Elastic" refers 
to the part of the response entering within the energy 
resolution width of the instrument. On IN5, this energy 
resolution was found to have a FWHM of AE ^120 /ieV. 
It can be seen from Fig. [13] that the elastic and inelas- 
tic intensities follow roughly the same Q-dependence as 
the static magnetic structure factor presented in sec- 
tion |III C| The intensity of the low-frequency response 
peaks at T ^ 70 K just above the macroscopic freez- 
ing temperature T g = 48 K (see inset of Fig. [l2| . The 
line-width remains small, of the order of V ~ 1 meV, 
independent of temperature. At the lowest tempera- 
tures T < 20 K the low-frequency inelastic intensity of 
magnetic origin essentially disappears - an observation 
backed up by our neutron spin-echo data shown below. 
This temperature corresponds roughly to the activation 
temperature of 25-30 K derived from the specific heat 
measurements shown in Sec. IIII A[ At the same time the 
elastic intensity increases smoothly below T « 80 K, 
indicating formation of a static magnetic ground-state. 
This is a rather standard signature of a dynamical mag- 
netic transition: a similar low-frequency response was 
also seen above the freezin g o f magnetic fluctuations in 
SrCr 9g Ga 12 _ 9g Oi9 (SCGOpl. 

We gain further insight into this phenomenon using 
a so-called elastic scan obtained on the backscattering 
spectrometer INI 6. We record the elastic scattering in- 
tensity I e \ during a temperature sweep over a Q-range 
of 0.25 A" 1 < Q < 2 A -1 . On IN16 the elastic energy 
resolution has a FWHM of AE w 1 /ieV. The elastic 
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FIG. 12: (top panel) The inelastic response of 

Y .5Ca .5BaCo 4 O7 at Q = 1.4 A -1 measured at 80 K 
on the IN5 spectrometer. The solid line shows a fit to a 
two-lorentzian model as described in the main text. The in- 
dividual broad and narrow lorentzian components are shown 
as dotted and dashed lines respectively. The instrumental 
resolution was modeled as a Gaussian of FWHM ^120 meV 
and is shown as the hatched area. 

(lower panel) Evolution of the inelastic response at 
Q — 1.4 A -1 as a function of temperature measured 
on IN5. The temperatures shown are 80 K (circles), 40 K 
(squares) and 10 K (triangles). The broad Lorentzian 
component is all that remains at low temperatures. The 
inset shows the temperature evolution of the elastic (squares) 
and inelastic (circles) integrated intensity at Q — 1.4 A -1 . 
For comparison, the elastic intensity measured at much 
greater resolution (FWHM ~ 1 /ieV) on IN16 is also shown 
(triangles). The scale is adjusted to highlight the relative 
difference with respect to data measured at the higher 
temperature range, in order to depict the development of the 
magnetic contribution. 
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FIG. 13: (Color online) The integrated "elastic" response (top 
panel) and integrated inelastic low-frequency (~ 1 meV) re- 
sponse (lower panel) of Yo.sCao.sBaC^O? derived from fit- 
ting spectroscopic data from IN5 to Eq. || as a function of 
temperature and wavevector transfer. 
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FIG. 14: (Color online) Upper panel - The intermediate scat- 
tering function on the correlation peak probed by NSE be- 
tween 2 K and 80 K. Lower panel - the fitted relaxation times 
(see text), an approximate Arrhenius behavior is indicated by 
the dotted line. 



and low-frequency inelastic intensities at Q = 1.4 A -1 
from both IN5 and IN16 are shown in the inset of the 
lower panel of Fig. [l2j The overall trend in the behavior 
of the elastic response on IN5 and IN 16 is very simi- 
lar, indicating only a weak dependence on the charac- 
teristic time-scale of the probe 51 , of the order of 10 9 Hz 
on IN16 and 10 11 Hz on IN5. Similarly, comparison of 



the temperature dependences in the inset of Fig. 12 with 
the total integrated magnetic moment measured on D7 
(which includes moments fluctuating at frequencies up to 
10 12 Hz) shown in the inset of Fig.pl shows that the spin- 
freezing behavior must be taking place over much longer 
timescales than those probed by the neutron measure- 
ments presented thus far. Accordingly it is interesting to 
compare these results with neutron spin echo (NSE) spec- 
troscopic measurements, which directly probe the time 
dependence over similar timescales. 

The data obtained with NSE show a continuous slow- 
ing down when cooling across the macroscopic freezing 
temperature T g . This has been observed in various spin 
glasses and is not too surprising given the different time- 



scales probed by the different techniques (/ ~ 10 9 Hz 
for NS E, H z < / < 10 3 Hz for DC or AC suscep- 
tibility).^^ At base temperature the response is truly 
elastic, S(Q,t)/S(Q,0) — 1 over the full time window 
probed. This is different with respect to the pyrochlore 
slab compounds SCGO and BSZCGO in which a fast ini- 
tial relaxation leads to S(Q,t)/S(Q,0) < 1, see Ref. [T71 
It also means that dynamics on the picosecond time-scale 
do not exist at all at base temperature which corrobo- 
rates the result hinted at by the apparent disappearance 
of the low-frequency Lorentzian component at low tem- 
peratures as seen from the IN5 data, suggesting the ab- 
sence of any magnetic fluctuations below fjw < 1 meV, 



see Figs. 12 and 13 



On increasing the temperature, the NSE data show a 
response that is much more stretched out in time than 
a single exponential relaxation function and gradually 
shifts down through the dynamic window of the instru- 
ment. This indicates that the response has contributions 
over many time scales simultaneously, making it difficult 
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to fit a single time scale (relaxation time) to the data, 
considering that the dynamic range of the instrument is 
short by comparison, covering less than three orders of 
magnitude in time. An attempt was made to quantify a 
characteristic time scale for the spin relaxation by fitting 
the NSE data to stretched exponential relaxation func- 
tion, 5(Q, i) oc exp — (t/r) . The resulting relaxation 

times are shown in the lower panel of Fig. [14] and ap- 
pear to follow an Arrhenius law: r = To • exp(A/T) with 
an energy barrier of A ~ 1580 K (140 meV), without 
any diverging behavior at T g . The phenomenological j3 
parameter is very small, j3 ~ 0.15 — 0.30, with a weak 
temperature dependence. 



IV. DISCUSSION 

The results shown above give some new in- 
sights into the frustrated magnetic properties of the 
Y .5Cao.5BaCo407 compound. The polarized neutron 
diffraction cross-section is in good quantitative agree- 
ment with the data reported earlier, 40 although the con- 
clusions we draw from our modeling of the data are some- 
what different. The most important new finding is the 
observation of an extended fast dynamic response that 
was missed up to now. This allows us to abandon the 
earlier hypothesis of magnetism based on the kagome lay- 
ers only, bringing our conclusions concerning the atomic 
structure and the magnetic short-range ordered struc- 
ture into line w ith other published work on undoped 
YBaCo 4 07l™ 

Our high-resolution neutron diffraction work reveals 
that the introduction of Ca into Yo.sCao.sBaCo^r re- 
sults in a disordered structure, with a splitting of the 01 
and 03 sites, and a consequent distribution of magnetic 
exchange pathways, depending on whether the C0O4 
tetrahedra are connected to Ca or Y ions. It is likely that 
this magnetic exchange disorder is responsible for the lack 
of long-range magnetic order in Yo.sCao.sBaCo^y. The 
symmetry of the crystal is trigonal and not hexagonal as 
was previously reported™ The introduction of atomic 
disorder and the fact that Ba 2+ is not under-bonded, 
may possibly be responsible for the lack of any signature 
of a trigonal-orthorhombic distortion in the magnetic sus- 
ceptibility of Yo.5Cao.5BaCo407 - although this has to be 
verified with low-temperature diffraction measurements. 

Sizable spin correlations develop below T ~ 100 K, 
leading to a well-developed short-range correlated but 
fluctuating state. The magnetic response remains fully 
dynamic down to T g « 50 K. At this temperature a par- 
tial freezing sets in as evidenced by bulk susceptibility 
and high-resolution inelastic neutron scattering measure- 
ments. Acknowledging the similar dynamic acceptance 
windows of D7 and IN11C (see section III]) we can con- 
clude that the diffraction response obtained on D 7 in the 
low-T < 20 K range is essentially elastic scattering prob- 
ing a static magnetic ground state. Accordingly the RMC 
spin configuration applies only to the frozen state. 



We have analyzed the experimental static structure 
factor S{Q) measured in the limited dynamic range of 
our cold neutron polarized diffraction measurement over 
the whole temperature range, from the low-T partially 
frozen state up to T = 300 K. Using the RMC method, 
we have examined various hypotheses concerning the lo- 
cation and presence of magnetic moments. Taking into 
account the deviation from ideal kagome geometry ex- 
hibited in the atomic positions of the (6c)-sites with ex- 
act and "kagome averaged" shell distances gives identi- 
cal results, suggesting that the bond length differences 
in the non-perfect kagome plane are not an important 
factor in determining the ground-state structure. The 
model restricting the magnetism to the kagome planes 
only results in a much poorer representation of the data 
than does the model with equal Co spins on all the Co 
sites. This indicates that there are no low-spin Co 3+ 
ions in Yo.5Cao.5BaCo407 as asserted by Schweika and 
co-workers. 40 Indeed it would be somewhat astonishing 
if low-spin Co ions existed within this crystal structure. 
The low-spin configuration of Co 3+ is not allowed in 
tetrahedral symmetry since it is not possible to have 
the two £2 electron spins antiparallel. A low-spin state 
would be stabilized if there were a strong site distor- 
tion which could perhaps lower one of the £2 orbitals 
to such an extent that it would be possible to over- 
come the Hund's rule coupling) 3 -^ From our diffraction 
data there is no evidence for such a strong distortion. 
Also, it was recently concluded from a combined exper- 
imental and theoretical x-ray absorption spectroscopy 
study of YBaCo407 that both the Co 3+ and Co 2+ are in 
the high-spin state in this material.^ Furthermore, the 
ideal stoichiometry of Yo sCao 5BaCo407 should be writ- 
ten (Y 3 +) .5(Ca 2 +)o.5(Ba 2 +)(Co 2 +) 2 .5(Co 3 +) 1 . 5 (0 2 -) 7 . 
Note that the Co 2+ to Co 3+ ratio is 2.5/1.5=1.67. In the 
P31c crystal structure, 3/4 of all Co ions form the kagome 
layered Co2 (6c)-site with the remaining 1/4 located on 
the Col (2a)-sites between the kagome layers. Thus the 
Co2 to Col ratio (=3) is higher than the Co 2+ to Co 3+ ra- 
tio (in YBaCo407 both ratios are equal to 3). This means 
that even if the Col (2a)-sites all have Co 3+ , the kagome 
layers must contain a distribution of Co 3+ and Co 2+ 
ions. Therefore we can assert that Y .5Cao.5BaCo407 is 
a mixed valence compound and the kagome layers must 
contain a distribution of mixed 5=3/2 and 5=2 spins. 
This mixed valence character adds to the disorder already 
present due to the Ca substitution and consequent ex- 
change pathway disorder. This reasoning highlights how 
the O-stoichiometry can influence the magnetic moment 
distribution on the Co sites and therefore play a critical 
role in determining the physical properties, in parallel to 
what was already concluded for YBaCo407 + 5 in Ref. [33l 
Oxygen stoichiometry and homogeneity may be the fac- 
tors explaining some of the differences seen in the bulk 
data of our sample when compared with that reported 
earlier W* 

The polarized neutron diffraction data further reveal 
that the total moment obtained over the restricted ex- 
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perimental dynamic range is considerably reduced with 
respect to the the expected full moment, in agreement 
with our observations of a wide dynamic range of mag- 
netic fluctuations using thermal neutron spectroscopy. 

The inelastic neutron measurements are consistent 
with fluctuating dynamics over a restricted spatial 
range with no signature of propagating fluctuations. 
The presence of propagating modes has been predicted 
for a classical kagome spin system^ - but only at 
low temperatures coinciding with the frozen regime of 
Yo.5Cao.5BaCo407. The dynamic response shows a com- 
plicated non-lorentzian line-shape. In the absence of a 
microscopic picture we have made an attempt to quantify 
the T-dependence by recognizing that at higher T the re- 
sponse is rather well described by a single lorentzian form 
with a width approaching r ~ 20 meV at T = 300 K. 
Below a temperature T ^ 100 K the intensity of the 
broad response flattens off, with a reduced line-width 
of T < 10 meV, however, at these temperatures a low- 
frequency response develops on top of the broad response. 
We obtain a picture in which a slowing down of dynamics 
precedes the eventual spin freezing and finally leads to the 
appearance of a fully elastic response. The appearance 
of this elastic part is seen in the high-resolution time-of- 
flight, backscattering and neutron spin echo data. The 
neutron spin echo results indicate a truly elastic response 
at base temperature, T = 2 K. We can conclude that at 
low temperatures the response of the spin system sepa- 
rates into two parts, one that freezes and another that 
remains dynamic. With the data available, however, it 
is not possible to identify an underlying mechanism for 
this on the atomic/molecular level. 

It is worth mentioning here that the magnetic 
properties of Yo.sCao.sBaCo^r bear some resem- 
blance to tho se of th e pyrochlore slab systems SCGO 
and BSZCGOEEEH but with some qualitative and 
quantitative differences. The dynamical response of 
Yo.5Cao.5BaCo407 shows rather clearly definable tem- 
perature dependent energy scales, in contrast to the 
cj/T-scaling and power law bulk susceptibility seen 
in the pyroch lore slabs and more recently in herbert- 
smithite! 10 * 56 ! It is also surprising that spin freezing in 
Y .5Cao.5BaCo407 happens at much higher tempera- 
ture in comparison with the overall energy scale. The 
role of disorder involved in the mechanisms leading to 
glassy ground state in the frustrated systems has at- 
tracted recent attentio n 1 57 * 58 1 and controlled studies on 



Y .5Cao.5BaCo407 might bring some more insight on this 
issue. 



V. CONCLUSIONS 

In this paper we have reported results and analysis 
of bulk and neutron diffraction and spectroscopic mea- 
surements made on the Yo.5Cao.5BaCo407 compound - 
proposed as a realization of classical kagome spin sys- 
tem. We have examined the magnetic signal of this com- 
pound over a broad dynamic range from below a fieV 
up to a few tens of meV. Our observation of a broad dy- 
namic response explains why the total moment carried by 
the Co-ions was earlier assumed to be too small, lead- 
ing to the conclusion that only the Co-ions residing on 
the kagome planes were magnetic. We have proposed a 
revised interpretation for the magnetic diffraction data 
obtained with cold neutrons, pointing to a situation in 
which all Co atoms are magnetic. In addition to the 
broad dynamic response with notable temperature de- 
pendence we have found conclusive evidence for a slow- 
ing down of part of the dynamic response, leading to a 
spin-freezing associated with the appearance of a truly 
elastic short-range correlated signal at the low tempera- 
ture limit. This spin-freezing is thought to come about 
due to the presence of exchange disorder between Co 2+ 
and Co 3+ ions in the vicinity of Y or Ca ions. We con- 
clude that Yo.5Cao.5BaCo407 is a highly frustrated anti- 
ferromagnet with an elastic response indicating a frozen 
component in the ground state. 
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